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Abstract 

In this paper we present a technique to measure the anisotropy parameter a suspension of 
microscopic particles. To estimate this parameter we use an integrating sphere as light collector, 
as well as a powder sample spherical particles of Ti02 with an average diameter a = 0.25 /im. 
The results for the sample of TiC >2 are consistent with Mie calculations. 
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I. INTRODUCTION 


The anisotropy parameter of a system of particles is defined as the average scattering 
angle, or as the first moment of the phase function of the system. The phase function of 
a system represents the probability density function of the scattered light. Therefore, the 
measurement of this parameter involves the implementation of an experimental arrange¬ 
ment that captures all the light scattered or as we will see, a ratio of the reflected and 
transmitted light. There are different methodologies to estimate this parameter, perhaps 
the best known is the method of scatterometer. 

The scatterometer is an optical instrument that allows us to estimate the scattered light 
as a function of scattering angle. However, one of its main drawbacks is that you need to 
implement a substantial infrastructure, plus a careful interpretation of the results. On the 
other hand, the method of integrating spheres allows us to implement a highly versatile 
experimental setup to measure this parameter. 

The integrating sphere method is commonly used to determine the optical properties of 
biological systems, its implementation usually requires the use of two integrating sphere. 
In the proposed settlement of this work will be done only using an integrating sphere and 
will be focused on the characterization of microparticles. 

The estimation of the anisotropy parameter through this method requires that the con¬ 
dition be fulfilled simple entertainment. As we will see, one way to do this is to estimate 
simultaneously the cross-section of absorption and extinction. 

The paper is organized as follows. In Section II it is present the theoretical framework 
that supports the experimental measurements. In Section III shows the experimental 
arrangements. In Section IV presents the materials and methods used in this study. In 
Section V presents the results and discussion of measurements. Finally, in Section VI 
presents the conclusions. 


II. THE ANISOTROPY PARAMETER 

The estimation of anisotropy (g) normally involves measuring the angular distribution 
of the scattered light. However, we will see that this parameter can be estimated using 
integrating spheres. 

For extremely diluted medium we can assume that the cells scatter from each other 
independently and there is no shading effects or multiple scattering. In this case, it is 
possible to match approximately the transmittance, reflectance and absorption properties 
of the medium with the cell anisotropy. Following is an analysis to justify this assertion. 
Suppose we have a medium with particles having a phase function P(6, <p). Must then be 
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satisfied that 


P{9,<f>)dn = 1 . 


( 1 ) 


' 47r 


If this particle distribution produces a pattern of isotropic scattering around (j), we have 


r2i r 


2n / P{9) sin 6d6 = 1. 


( 2 ) 


With the change of variable x = cos 9 the integral is simplified and we can write 


p(x)dx = 1, 


( 3 ) 




where the function p(x) = 27rP(cos 1 x). 

The anisotropy parameter is defined as the first moment of the phase function. That is, 


g= xp(x)dx. 


( 4 ) 




Since x = cos 9, g represents the average cosine of the scattering angle (as already men¬ 
tioned). 

The fraction of scattered light that is scattered forward is given by, 


/ = / p(x)dx. 


( 5 ) 


Similarly, the fraction of scattered light exiting backwards is 


b= p(x)dx. 




and should to comply that / + b = 1. 
We write now, 


g = gg + gb 


where 


and 


gf= xp(x)dx 


gb = L xp{x)dx. 


( 6 ) 


( 7 ) 

( 8 ) 
( 9 ) 


We consider first the case of large particles. In this case, the spread is primarily at small 
angles, either forward or backward. For these angles x = cos 9 ~ ±1 , so gf ~ / and 
g b ~ -b. 
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With this approach we need 
Using the condition, 


f-b = g. 


f + b = l, 


( 10 ) 


( 11 ) 


and solving for / and g , we have that the fraction of light scattered forward can be written 
as 


/ = 


i + g 
2 ’ 


while the fraction of light scattered backwards is, 

6 = 1 ^. 


( 12 ) 


(13) 


We can see that a value of g = 1 (all light scattered forward) results in / = 1 and 6 = 0. 
Furthermore, with g = —1 (all the light scattered backwards) we obtain / = 0, and 6 = 1. 
This is consistent with expectations. 

Furthermore, for small particles it has isotropic scattering. In this case, p(x) = p(—x), 


g = 0, and f = b = 1/2. We see that this result also agrees with the approximation (12) 


and (13). This suggests that those who were motivated by an approximation for large 


particles, may be reasonable approximations for particles of other sizes. This suggests 
that these approaches, who were motivated by an approximation for large particles, may 
be reasonable to particles of other sizes. 

It is clear that the validity of these expressions for intermediate values of g depends on the 
phase function, but it is interesting to explore if the phase function of Henyey-Greenstein. 
Figure [I] shows the behavior of the values of / and b calculated using the expressions (12) 


and (13) (solid line) and calculated based on the phase function of Henyey-Greenstein 
(dotted line). As expected, the curves coincide in the limit g — — 1 y <7 = 1, as well as 
the center point g = 0. There are only small deviations for other values of g. This gives 


some confidence in the approach represented by the expressions (12) and (13), which also 
implies that g = f — b. 

This shows the possibility of estimating the anisotropy factor g in terms of the properties 
of reflectance and transmittance of the diluted medium. That is, without requiring angle 
scattering measurements. 

We know that, by definition, absorption (A), reflectance ( R) and transmittance (T) satisfy 


A + TZ + T = 1 

Another way to write this relationship is, 

r n 


1 -A 1 -A 


= 1 . 


(14) 


(15) 
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FIG. 1. Comparison between the values of / and b with approximation x ~ ± 1 (straight line) 
and the phase function of Henyey-Greenstein. 


The reason for doing this is to make clear that this is equivalent to the condition f + b 
so we have to 


/ 


T 

1-A’ 



n 

^~A' 


= 1 , 

(16) 


Hence, it is possible to determine the value of g through measurements of transmittance 
and reflectance coefficients. Thus from equations (10) and (16), 


T-n 
r + iz' 


(U) 


Note that for this we are assuming that media is sufficiently diluted as to have only single 
scattering, plus you can collect all the transmitted and reflected light. 

We can say that since we have 7Z, T and A, with the same arrangement we can estimate 
the parameter of anisotropy of the cells by dilutions of the medium. Although in principle 
we would have spurious contributions in reflectance due to air-interface cell plastic we can 
correct it by using a reference cell and removing the specular component measurements. 
The arrangement allows us to estimate the absorption (^4) from the reflectance measure¬ 
ments ( 1Z ) and transmittance (T) of the medium, as mentioned above. 

Once we have estimated absorption, we used the relationship between this and the ab¬ 
sorption cross section C a of the particles, 


A = exp(-pC a L), (18) 

where A is the ratio of the the reference power P 0 (measured in the absence of particles) 
and the reduced power P' n (with the particles present), p is the density of particles, and 
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L is the thickness of the cell [see equation (4) from Ref. [5j]. Solving for C a , we have, 


C a — — — ln(^4), (19) 

7] 

where rj = pL represents the number of particles per projected area. 

As we will see in the next section the estimation of this parameter allows us to determine 
the anisotropy of the particles adequately. 


III. EXPERIMENTAL SETUP 


The diagram in Figure [2] shows the way in which we estimate the amount of scattered 
light using an integrating sphere as a collector of light, a quartz cell with the particle 
suspension and a photodetector for measurement. It is noteworthy that for particles with 
diameters in the micron range most of the scattering is forward, so we can consider that 
the light returning from the field to the sample does not influence the measurements, 
which simplifies the estimation of the amount reflected or transmitted by the suspension 
of light particles. 

The arrangement consists of an integrating sphere Ocean Optics (Model ISP-50-8-R-GT), 
a NewPort (Model 2051-FS-M) photodector, a laser HeNe (12 mW, random polarization) 
a quartz cell with 1 cm thick, 5 permanent mirrors and 2 removable. The right arm of the 
array illuminates the particles in transmittance and the left arm in reflectance. Note that 
disabling the removable mirror 1, the beam illuminates the sample directly to an angle 
8 °, the integrating sphere has an aperture placed in this position to pass up the specular 
component. On the other hand, the beam 3 also illuminates the sample at an angle 8° 
transmittance but in this case we also let the specular component leave the sphere. The 
integrating sphere employed (Labsphere 4P-GPS-053-SL) has a wall reflection coefficient 
m = 0.99, an internal diameter of 13.46cm, two 2.54cm diameter ports, and a fiber-coupled 
exit port in which we connected the fiber optics photodetector (?). 

The photodetector will collect a detected power Pd when illuminate with the beam 1. 
Then if we illuminate with the beam 2 we will get a power P r . However, if we illuminate 
with the beam 3 we will get a power P t . From Ref. [5] [see equations (8) y (9)], we have 
that the reflectance and transmittance can be written as 


n = m— and 
Pd 


q- Pt 
T = m—. 
id 


( 20 ) 


From the equations ([17]) and (20), we have that the anisotropy parameter can be written 


as, 


9 = 


Pr-Pt 
Pr + Pt 


( 21 ) 
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FIG. 2. Scheme of the experimental setup to estimate the parameter of anisotropy of micro¬ 
particles. Figure (a) shows the configuration to determine the proportion of light scattered 
forward /. Figure (b) shows the configuration to determine the proportion of light scattered 
backward b. 

The method was applied to a sample of powdered titanium dioxide (TiC^). The main 
advantage of working with this kind of particles is that we can characterize optically with 
Mie theory, which allow us to check the consistency of the experimental measurements. 
Another way to check the consistency of the experimental measurements is through the 
estimation of the absorption and the extinction of the sample. In order to do this we use 
the methodology presented in Ref. [5J. The arrangement of the Fig. [3f^a) allows us to 
estimate the extinction cross section C e through the expression 

P m = exp (~pC e L), (22) 

where P D is the reference power (without particles), P m the power collected by the pho¬ 
todetector (with a suspension of particles), p is the density of particles, and L is the 
thickness of the cell. On the other side, the Fig. |3](b) allows us to estimate the absorbtion 
cross section C a through the expression 


P' m = P 0 exp(-pC a L ) 


(23) 
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FIG. 3. Schematic diagrams of the arrangements used for the characterization of the photosyn¬ 
thetic cells. Setup for the measurement of (a) extinction and (b) absorption [SJ. 


where P. [ ' m the power collected by the photodetector (with a suspension of particles). We 
should note that although the equations (22) and (J23| are very similar estimation of the 
values of P m and P' m are very different. For a deeper description of the estimation of these 
parameters see Ref. |,5j. 

We did dilutions of the sample to verify the linear dependence between the number of 
particles per area projected rj and the cross section of absorption C a , which is true only if 
the condition of independent and single scattering are valid. Therefore, we consider that 
the measurement region for estimating valid absorption cross sections and extinction of 
particles is in the range where the measurements follow a linear behavior as a function of 
the dilutions, see Fig. [IJ 


IV. RESULTS AND DISCUSSION 

In this section we estimate the anisotropy of the sample of Ti02 and we compare this 
value with that obtained with the Mie theory. The sample of TiC >2 that we consider is 
composed of particles of an average diameter 0.25 /jrn, an index of refraction of particle 
rir = 2.093 immersed in water. In this case the Mie theory leads us to obtain an anisotropy 
parameter g = 0.59 and a extinction cross section of C e = 0.068 /irrr to A = 632.08 nm. 
We will start our measurements ensuring that the condition of single scattering is satisfied. 
So it was necessary to determine the absorption and extinction spectrum of the sample. 
Experimental measurements of absorption and extinction of Ti02 are shown in Figs. 
[4j The Fig. |4|a) was obtained with the arrangement of the integrating sphere Fig. [2] 
and Figs. [I] (b) and [I] (c) the absorption and extinction arrangements shown in Figure 
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FIG. 4. Parametro de anisotropfa de polvo de TiC >2 de partfculas de 0.25 pm de diametro para 
diferentes diluciones. 


??, respectively. In this case four sample dilutions were made, dilution 0 represents the 
original medium, dilution means 1 2 diluted in parts and so on. We can see that to 
dilution 0 the media was highly concentrated so that we started from the dilution 1 -see 
Figs. 3 b ) and [4] (c)- and ended up in the third dilution because the environment was 
very diluted. The dotted red line represents a least squares fit of the selected points. 

The absorption cross section C a of Ti0 2 can be estimated through a least squares fit to 
the points shown in Fig. |4|a) and equations([l9|), 

> = 1.16 x 1(T 3 pm^ 2 . (24) 

On the other hand, from Fig. [4](b) we have that the absorption cross section obtained 
with absorption arrangement is Ref. [5], 

C¥ ] = 2.95 x 10 -2 pm" 2 . (25) 

where the number of particles per projected area 


rj = 96.6 pm, 


(26) 


was obtained from the experimental measurements of extinction [Fig. |4](c)] and the ex¬ 
tinction cross section C e obtained with Mie theory. 

We can see that there is a significant difference (one order of magnitude) in the estima¬ 
tion of this parameter, if we use the integrating sphere arrangement or arrangement of 
absorption. This is due to that the arrangement of the sphere captures the light coming 
forward and back, while the arrangement of absorption captures light only going forward, 
which leads us to have an overestimation of C a , as we can see from the values shown in 


equations (23) and (24). However, for more absorbent or larger particles the difference is 


not so marked, because there will be less light scattered backwards, especially for spher¬ 
ical particles, however, for particles with other geometries the situation may change and 
you need to take into account the limits of the array of absorption. 
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FIG. 5. Anisotropy of TiO 2 as a function of the dilution of the media. 

The anisotropy of the particles of TiC >2 was obtained from the values shown in Figure [5j 
The vertical axis shows the measurements obtained with the arrangement of the sphere 
(see Figure |2| and equation (??), the horizontal axis the dilutions of the media. The dilu¬ 
tion 0 corresponds to the media undiluted where 7] = 99.6 [i m -2 , each dilution represents 
the third of the above, so that dilution for 2 we have rj/9 = 11.06 /im~ 2 . 

We can see that the variations in the anisotropy parameter is related to the cases where 
we have many or too few particles. Remember that the anisotropy parameter is valid 
only when we have single sacttering, like C a , so that the region where the latter is valid 
it will also g. Thus, in Fig [5j we will consider as valid those points where the value of 
g is constant or take their greatest value this case corresponds to dilution 2 , so that the 
anisotropy of Ti 02 is, 

g = 0.615. (27) 

As we have seen it is very important to ensure that the condition of single scattering is 
fulfilled. One way to know this is through the product of 77 x C s , which gives us an estimate 
of the ’’area” occupied by the particles in the illuminated region. For values of rjC s near 1 
or higher imply that multiple scattering events occur, while for very small values rjC s <C 1 
will have single scattering. In the case, we have rj/9 x C s = (99.6/9) x (0.06684) = 0.73, 


TABLE I. Propiedades opticas del TiC> 2 . 



C e [fjm 2 ] 

C a [/nn 2 ] 

g H 

Experimental 

- 

1.16 x 10" 3 

0.615 

Calculado 

0.068992 

- 

0.59945 
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which tells us that dilution 2 is a good approximation to single scattering. 


V. FINAL REMARKS 

The results show that for the particles of TiO -2 the experimental and theoretical results 
are consistent. However, the borosilicate glass powder has inconsistencies between these 
results. This may be because in principle we would have spurious contributions in re¬ 
flectance due to the air-plastic cells, we could be corrected using a reference cell. 
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